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It is known that baryon number inhomogeneities may arise as a consequence of 
- electroweak baryogenesis. Their geometry, size, and amplitude depend on the pa- 

' rameters that characterize the baryogenesis mechanism, as well as on those that 

■ determine the phase transition dynamics. We investigate this parametric depen- 
^ . dance. We show that in the case of the minimal supersymmetric standard model, 
Q ! the geometry of the inhomogeneities most probably consists of spherical regions of 

Q I high density surrounded by low-density walls, in contrast to the case of the minimal 

psj ' standard model. In this supersymmetric extension we find that density contrasts of 

■ up to a factor of 100 may arise. This amplitude increases for higher values of the 
. latent heat or lower values of the bubble wall tension, and can be significantly larger 

^ , in different extensions of the standard model. Such inhomogeneities may thus affect 

I the dynamics of the subsequent quark-hadron phase transition. 

cn ■ 

a^ : 

O . I. INTRODUCTION 

■ 

o : 

^ I The generation of inhomogeneities in the barvpn number density is a possible outcome of 

P^l the QCD and the electroweak phase transitions [l|, Q, 0, 0| • These may arise as a consequence 
p • of the perturbations produced by the walls of expanding bubbles in the surrounding plasma. 
QJ . In the QCD case, baryons, which are much heavier in the hadron phase than in the deconfined 
quark phase, are pushed away by the walls of expanding bubbles as the phase transition 
, develops. Then, as the volume occupied by the hadron phase grows at the expense of the 
^ ■ quark phase, baryons are driven into small regions of space P. Hence, it is believed that 
the general geometry of the QCD inhomogeneities is that of localized clumps of high density 
surrounded by voids of low density. 

In the electroweak phase transition, baryon number inhomogeneities can arise as a con- 
sequence of the generation of the baryon asymmetry of the Universe (BAU), whether the 
latter is produced through the standard mechanism of electroweak baryogenesis |3] or in 
electroweak cosmic strings yj. In general, the electroweak phase transition provides the 
three Sakharov conditions p for baryogenesis, namely, baryon number violation, C and 
CP violation, and a departure from thermal equilibrium. These conditions must combine 
in such a way as to produce a net difference between number densities of baryons and an- 
tibaryons. Interestingly, the idea of electroweak baryogenesis requires physics beyond the 
minimal standard model (SM) in order to give a quantitatively successful result. The basic 
idea is the following (for reviews on electroweak baryogenesis see @|). 

The standard mechanism assumes a first-order phase transition in order to achieve the 



JS 



non-equilibrium condition. As bubbles expand, a chiral flux is injected in front of their 
walls due to CP violating interactions of the bubble walls with the particles of the plasma. 
Thus, an asymmetry between left handed quarks and their antiparticles is generated near 
the interfaces. This asymmetry biases the baryon number violating sphaleron processes in 
the symmetric phase. The resulting baryon asymmetry is caught by the walls and enter 
the bubbles. In order to avoid the washout of the generated BAU when equilibrium is 
restored, the sphaleron processes must be suppressed in the broken symmetry phase. This 
requirement imposes a condition on the value of the Higgs fleld (j) inside the bubbles 0, 

0„ (T) /T > 1 . (1) 

Here, 0^ is the vacuum expectation value (VEV) of at temperature T, which corresponds 
to the global minimum of the free energy. It plays the role of an order parameter, and 
the condition ((T)) states that the phase transition must be strongly flrst-order. This usually 
constrains the parameters of the theory. 

The generated BAU has a strong dependence on the velocity of bubble walls. On 
one hand, if the velocity is too large, the left-handed density perturbation will pass so 
quickly through a given point in space that sphaleron processes will not have enough time 
to produce baryons. Thus, the resulting baryon number density hb will be small in this 
case. On the other hand, for very small velocities thermal equilibrium will be restored and 
the baryon asymmetry will be erased by sphalerons. Hence, the baryon production will be 
small again. Therefore, the generated baryon number has a maximum for a certain wall 
velocity = Vp^ak, which can be estimated by comparing the baryon number violation time 
scale with the time of passage of the chiral asymmetry Such estimates give a small 

value fpeak ~ 10~^, which is conflrmed by numerical calculations 0, U . 

This dependance oi ns on the bubble wall velocity is important for electroweak baryo- 
genesis, since it establishes a preference for models in which the velocity is close to fpeak- 
On the other harid, the velocity of bubble expansion is not constant throughout the phase 
transition llM ll^ ■ The variation of causes a variation of the local baryon density 

left behind by the walls along the bubble radius. It is thus evident that generating the ob- 
served BAU in the electroweak phase transition entails the formation of inhomogeneities in 
the baryon number density. These inhomogeneities are spherically symmetric and centered 
at the nucleation points. 

A general feature of flrst-order phase transitions is the slow-down of bubble expansion 
due to the release of latent heat Jl3']. In the case of the electroweak phase transition, it is 
known that the velocity of bubble walls may decrease from an initial value Vi ~ 10~^ — 10~^ 
to a minimum velocity Vm ~ 10^^ — 10~^, but the variation can be larger, depending on 
the model 0, 3|. The maximum variation of the baryon density ns as v^, decreases 
from Vi to Vm determines the amplitude of the baryon inhomogeneities generated inside each 
bubble. The exact proflle nsir) along the radial direction and the general geometry of 
the inhomogeneities depend on the evolution of For instance, if the initial velocity 

Vi is close to the peak velocity fpeak, then clumps of high baryon density will be formed, 
surrounded by voids of low density. On the contrary, if Vm — "Wpeak, walls of high density will 
surround voids of low density. The size of the inhomogeneities is roughly given by the mean 
separation between the centers of bubbles at the end of the phase transition. 

The survival of the inhomogeneities at later epochs depends on the amplitude of the 
fluctuations as well as on the separation of the centers of fluctuations, which must be larger 
than the diffusion length of baryons for that fluctuation amplitude. In general, the inho- 
mogeneities will have different cosmological consequences, depending on the epoch in which 



the phase transition occurs. Baryon inhomogeneities generated in the QCD phase transition 
have been intensively investigated due to their influence on Big Bang Nucleosynthesis (BBN) 
P, 0, 0] . On the contrary, baryon inhomogeneities generated at the electroweak scale are 
not expected to survive long enough to affect BBN [Jl^- 

However, electroweak inhomogeneities may survive until the QCD scale 0, 12 1. Then, 
the dynamics of the quark-hadron phase transition may be affected by their presence 
0]. This is because the critical temperature Tc of the quark-hadron phase transition depends 
on the chemical potential, and therefore on the number density, being lower in regions 



with higher baryon density |j2CI|. I21|. Hence, the process of bubble nucleation is delayed in 
such regions. The effect of pre-existing baryon inhomogeneities on the quark-hadron phase 
transition depends on the profile of the inhomogeneities. 

Assume for instance that the geometry is that of small regions with low baryon density 
in a background of higher density. Then, since the critical temperature is reached first 
in the low-density regions, the inhomogeneities may act as impurity sites where hadron- 
phase bubbles can nucleate before the supercooling required for spontaneous nucleation is 
reached flS*]. Thus, the QCD phase transition may proceed by inhomogeneous nucleation 
at Tc, rather than by supercooling and homogeneous nucleation with a rate F per unit 
time and volume. Therefore, the size scale of the inhomogeneities determines the mean 
separation between centers of nucleation. This distance, in turn, sets the scale of the QCD 
inhomogeneities. 

If, on the contrary, the inhomogeneities have the form of small lumps of high baryon den- 
sity surrounded by a background of lower density, then the phase transition will be delayed 
inside the lumps. Since the Universe quickly reheats close to Tc, the high density regions stay 
in the quark-gluon plasma phase, while the hadron-phase bubbles grow outside. This also 
affects the QCD generation of inhomogeneities, since baryon number tends to accumulate in 
the deconfined quark phase, as explained above. As a consequence, inhomogeneities gener- 
ated at the electroweak scale may be amplified at the QCD phase transition, and therefore 
indirectly affect primordial nucleosynthesis [l6| . 

In this paper we study the formation of baryon density inhomogeneities in a baryon- 
generating electroweak phase transition. This issue has been previously investigated by 
Heckler j2| for the minimal standard model. However, as we already mentioned, electroweak 
baryogenesis is possible only in extensions of the SM. The function ub (vw), as well as the 
dynamics of the phase transition, can be quite different in such extensions. As a consequence, 
the kind of inhomogeneities that form may be quantitatively as well as qualitatively different. 
In particular, in Ref. [3| the wall velocity was supposed to lie well on the right of tWk, and 



a dependance ub oc was assumed. In contrast, recent calculations [10|, |ll|, show 
that the initial velocity is likely to be around Wpeak, at least in the case of the minimal 
super symmetric standard model (MSSM). In such a case, the velocity slow-down may occur 
completely on the left of the peak |U] and, as we shall see, yield the opposite behavior, 
Ub oc f^. 

In sec. |H]we analyze the dependance of ub on Vw in the whole range of (non-relativistic) 
wall velocities. In sec. II I II we study the evolution of the wall velocity using analytic approxi- 
mations which allow to determine the parametric dependance of the amplitude and shape of 
the inhomogeneities. Unfortunately, such approximations give only a qualitative picture, so 
we also make a more numerical investigation of the electroweak phase transition: In section 
II VI we derive appropriate equations for the evolution of temperature and fraction of volume 



occupied by the broken-symmetry phase. We integrate this set of equations numerically and 



calculate the baryon number density. In section |V| we present the results of this calculation. 
We consider different values of the parameters, which include the case of the MSSM. Finally, 
in section IVII we discuss the possibility that the electroweak baryon inhomogeneities have 
an effect on the QCD phase transition. We argue that this is possible in some extensions of 
the SM. Our conclusions are summarized in section IVIII 



II. DEPENDANCE OF ELECTROWEAK BARYOGENESIS ON THE WALL 

VELOCITY 

The baryon asymmetry originates from the chiral quark asymmetry in the symmetric 
phase that is generated as a consequence of the CP violating currents coming from the 
wall. The density of left-handed particles in front of the bubble walls can be calculated by 
considering a set of coupled diffusion equations of the form (see e.g. 0,01) 

= A V'n, - J2 Ki^a^f^ + (2) 

for all particle species i that have CP violating interactions with the bubble wall. Here, Di 
are the diffusion constants, are the rates of reactions that change the number of particle 
i by an amount A^j (the corresponding rates per unit volume are T^F^/G), ki are statistical 
factors: ki = Qi for fermions, ki = 2gi for bosons, where Qi is the number of degrees of 
freedom of species i, and 7^ are the rates of generation of axial numbers at the wall, which 
depend on the Higgs profile. 

In general, the diffusion equations (j2I) must be solved numerically. However, several ap- 



proximations can be made in order to simplify the problem [10, |11|, |2J, |2^. For instance, 
the diffusion constants for all quarks can be assumed to be nearly the same, and some of 
the rates F^ (e.g., those corresponding to strong sphalerons and trilinear Higgs-quark inter- 
actions) can be assumed to be fast enough to be in thermal equilibrium. As a consequence, 
several of the particle densities can be eliminated algebraically, leaving only a few equations 
to deal with. In order to avoid entering into the details of any specific model, we will assume 
that all the Ui can be written in terms of just one of them, so only one diffusion equation 
needs to be solved (see e.g. [IS])- This will allow us to find an analytical solution for the 
baryon density, which can be easily inserted in the numerical calculation of the evolution of 
the phase transition. 

Therefore, we can write a diffusion equation for the density of left-handed fermions with 
effective values of D, F, and 7, 

Dnl + vWl - Fn^ + 7 = 0, (3) 

where we have made the usual assumption of a planar wall which moves non-relativistically 
to the right, so 7^ and rii depend only on 2; — Vy^t. The diffusion constant D for the chiral 
asymmetry depends on the particle spectrum of the specific model. In supersymmetric 
extensions of the SM, the main contributions to D come from the Higgs sector and give 
D ~ 100T~^ ji^]. The effective axial quark number relaxation rate F receives contributions 
from the helicity flip due to the t op q uark mass and from Higgs self-interactions. In the wall 
frame, it can be approximated by Ifl U, l^] F (2) = TO {—z), where 9 is Heaviside's function 



and f ~ 10 ^T. The weak sphaleron processes are very slow, so they can be disregarded 
for relatively large wall velocities. However, we will consider a range of velocities that may 



include very small values of Vw Therefore, we should take into account that the effective 
rate F in the symmetric phase is not zero but of the order of the weak sphaleron rate. 
Consequently, we will assume that F = F in the broken-symmetry phase, and F = aV^g 
in the symmetric phase, where a is a numerical constant that depends on the degrees of 
freedom contributing to Eq. (jH)), and F^^ ~ 10~^T ji^. 

Thus, Eq. (jni) splits into two identical equations for 2; < and z > The solution is of 
the form 

riL = Ae-'^+" + Be-^-' (4) 
- [D {io+ - io-)]-^ r [e~"-(^-^') - e-"+(^-"')l 7 (/) dz', 







where u±, A and B are different on each side of the wall. The parameters uj± are given by 

u± = v^/2D ± ^Jiv^/2Df + T/D, (5) 

and the constants A and B are determined by the boundary conditions rii (±cxo) = and 
continuity of and n'^ aX z = 0. Their general expressions are rather cumbersome and can 
be found in the appendix. If the CP violating source 7 is localized on the wall, then outside 
the wall riL (z) has a simple exponential dependence [notice that, since Ci;+ > and < 0, 
only one exponential function survives on each side of the wall in Eq. 

Since baryon number violation takes place in the symmetric phase, we are only interested 
in the solution for z > 0, 

riL (z) = Ae--^+^ (6) 

where Us+ is given by Eq. (0), with F = aT^s, and the constant A is an integral of 7 (z) 
times a combination of exponentials. The source 7 is proportional to v^, since it is given 
by the flux of particles that reflect from the moving wall '26]. Thus, A is the product of 
times a factor that depends on the details of CP violation at the bubble wall. Although this 
factor can be velocity-dependent, in general it is not very sensitive to Vw In the appendix 
we check this assertion using an approximation for the source 7 (z) . According to these 
considerations, we will assume that A oc Vw Since we are interested in density contrasts 
caused by a variation of the wall velocity, we only need to determine the baryon density ub 
up to a constant factor independent of (we assume that the CP violation is enough to 
generate the observed BAU). 

The baryon number density also satisfles a diffusion equation, similar to Eq. 0. In this 
case, only sphaleron processes are relevant, so we have 

Dgu'^ + v^n'j^ - 3T^sO (z) {ul {z) + feng) = 0. (7) 

Here, Dq ~ 6T~^ is the diffusion coefficient for quarks, and the numerical factor h depends on 
the particle spectrum. The terni proportional to Ub accounts for baryon number relaxation 
when the wall velocity is small (ifll. The chiral density ni acts as a source term in Eq. 
dZj) and is given by Eq. (jH)). 

The solution to Eq. (jZj) for z < is a constant, ub = ub (0). For z > 0, has the form 
of Eq. (jH), with 7(z) = 3F^s'n,i(z) and 



UB± = v^/2Dg ± ^J {v^/2Dgf + SbT^s/D^. (8) 



The conditions of continuity at z = and vanishing baryon number aX z = +00 (i.e., far in 
the symmetric phase) give the value of the constant baryon density in the broken symmetry 
phase, 

nn = / nr. {z) e^^-^dz. (9) 

The parameter ujb- is often approximated by ojb- — —'^hVy^s/vw, which is correct for 3> 
ADqVyjs, i.e., for Vyj > 10^^. Eq. Q can be easily integrated using Eq. 



riB 



(10) 



where C is a constant. 

The analytic approximation (|TO|l for the dependance of the BAU on the wall velocity is in 
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2^ 
of a 



reasonable agreement with other approximations and numerical calculations [1 
Two quantitative differences with some works arise due to our inclusion in Eq. (jHI 
non-vanishing rate F ~ F^g in the symmetric phase, and the inclusion of the diffusion 
coefficient for quarks in Eq. ((7j). In Fig. [T]we plot {vw) for a = 36 = 1. We see that 
the baryon density has a peak, and tends to zero for small or large wall velocities^. For 
Dg ^ D, the peak is close to = a/DT\^. This can be seen either numerically or by 
making approximations of Eq. pO|) for the different possible ranges of v^. Furthermore, it 
can be easily seen that Eq. ()10|1 has the expected behavior on each side of the peak, namely, 
ub oc Vw for <^ V-DF^s, and ub oc t;^^ for v.^, > V-DF^s- 
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FIG. 1: Baryon number density as a function of the wall velocity 



As we shall see, the wall velocity decreases during the phase transition as a consequence 
of latent heat release. Therefore, inside a given bubble, the local value of the baryon number 
density at a distance r from the center of nucleation is given by ^.^(r) = [^^«)(^)], where t 
is the moment at which the bubble wall has passed through r, i.e., t is given by R{t) = r. 
The total BAU is the volume average B = / [fu, (t)] df, where / {t) is the fraction of 
volume occupied by bubbles. We remark that in this work we are assuming that B gives 
the observed baryon abundance, which depends on the value of the constant C in Eq. (jlOj) . 
We are interested in the relative variation of the baryon number density left behind by the 



Notice however that this approximation breaks down for v^, close to 1. 



bubble walls, which does not depend on this constant, and can be computed from Eq. ()10|) 
once the time dependance of the wall velocity is known. Still, before closing this section we 
would like to comment on the effects of the velocity decrease on the total amount of baryons 
B. 

If the initial velocity is Vi ^ fpeak ~ VdT\^, a decrease of produces an enhancement 
of the BAU On the contrary, for Vi -C fpeak, a velocity decrease will cause a suppression 
of the generated baryon asymmetry 14]. It is interesting to notice that in the case of 
an enhancement the effect can be quite large; however, in the case of a suppression the 
electroweak baryogenesis scenario will not be harmed significantly. Indeed, the volume 
spanned by a bubble wall while the velocity is closer to the initial value Vi is in general of 
the same order of that spanned with closer to the minimum velocity Vm- Assuming that 
hb takes very different values hb (vi) and ns (vm) in each of these volumes, then the total 
baryon number density will be of the order of the largest one. Even in models in which 
baryon densities with opposite sign are formed in different regions, the suppression to the 
total BAU should be of order one. Otherwise it would require a fine tuning between the 
dynamics of the phase transition and the baryogenesis mechanism. 



III. PHASE TRANSITION AND BARYON INHOMOGENEITIES 

The electroweak phase transition takes place in the radiation-dominated era, at temper- 
ature Tc ~ lOOGeV, when the expansion rate is He = l/2tc and the age of the Universe is 
given by 

= iMjTl (11) 

Here, Mp is the Planck Mass, and the parameter ^ is related to the light degrees of freedom 
of the plasma, (7*, by ^ = ^/9^J^2tt^^ (at the electroweak scale, ~ 100 and ^ ~ 1/34). 
The characteristics of the transition are determined by the free energy density difference 
V (0, T) between a value 7^ of the Higgs field VEV (corresponding to the broken- 
symmetry phase) and = (symmetric phase). In a first-order phase transition, V" (0, T) 
has a non-zero minimum 0^ {T) coexisting with the minimum = and separated from it 
by a barrier. Then, the quantity V (T) = V (0m (T) , T) gives the pressure difference between 
the two phases. In the simplest case, four parameters characterize the behavior of V (0, T), 
namely, the critical temperature Tc at which V (Tc) = 0, the value of the order parameter at 
Tc, 0m (Tc) /Tc, the latent heat (i.e., the energy density discontinuity at T = Tc), given by 

L = Tc {dV/dT)j,^ , (12) 

and the tension of the bubble wall, 




where (r) is the wall profile at T = Tc. 

An important parameter, which is not directly related to the free energy but is relevant 
for the dynamics of the phase transition is the friction coefficient 77. In a first-order phase 
transition the system supercools, so bubbles nucleate at a temperature T/v < Tc. At T < T/v 
bubbles expand with a velocity that is determined by the pressure difference between the 



broken-symmetry phase inside the bubbles and the symmetric phase that surrounds them, 



and by the friction of the wall with the plasma [8|, |22, |27|, |2 



(T) = -V (T) /r/. (14) 

The electroweak wall velocity is subsonic, i.e. < c^, where = 1/a/3 is the velocity 
of sound in the relativistic plasma. Therefore the bubble expands as a deflagration. The 
deflagration front is the bubble wall, where the value of changes from to 0^, and the 
energy ()12|1 is released. A shock wave is formed in front of the wall, which is preceded by a 
supersonic shock front [1^ . The energy released in the deflagration process is transmitted to 
the fluid and distributed between the two fronts. Most of it is closer to the phase transition 
front than to the shock front 0, H^, so a temperature proflle arises. However, if -C Cg 
(which is the case of the electroweak bubbles), part of this heat is carried away far in front 
of the wall and can influence other bubbles. If v^j is small enough, one can assume that the 
released latent heat has enough time to distribute uniformly and equilibrate the temperature 
everywhere |3|. We will thus assume for simplicity that the only effect of the shock wave is 
to cause a uniform reheating of the plasma. This will be a good approximation if the time 
it takes the latent heat to get distributed throughout space (which depends on the bubble 
separation) is much shorter than the time scales involved in the development of the phase 
transition. We will examine the consistency of this approximation in section |3 

The phase transition thus develops essentially in two steps (see e.g. [3]). The flrst stage 
is characterized by a quick increase of temperature and a decrease of f^, as the latent heat 
reheats the plasma. After a certain time Sti the plasma has reheated to a temperature that is 
close to the critical one, and the system enters a phase-equilibrium or slow-combustion stage 
ll^]. During this second stage the latent heat release is compensated by the expansion of 
the Universe, and the temperature remains almost constant. Since T ~ Tc, (T) ~ 0, so 
bubble expansion slows down signiflcantly, and bubble nucleation effectively stops. Hence, 
this stage lasts for a longer time 6t2. 

The time 5tr during which bubbles nucleate is much less than 5ti and St2 0]. This 
is due to the quick variation of the nucleation rate F (T) with temperature. Therefore, 
it is a good approximation to assume that all the bubbles nucleate at the beginning of 
the phase transition (this is conflrmed by our numerical calculation). Hence, the number 
of bubbles remains flxed during the phase transition, and their number density can only 
change due to the dilution caused by the expansion of the Universe. However, the time scale 
of the electroweak phase transition is much less than the age of the Universe, so this effect 
is negligible. This means that when the phase transition completes, all the bubbles have 
roughly the same size, which is ~ nb(Tj 



N 



-1/3 



The size scale of the inhomogeneities is given by the size of the bubbles, and hence by the 
distance between centers of nucleation. Their shape depends on the distances travelled by 
the bubble walls during the times 6ti and 6t2, and on the baryon- number density generated 
during each stage. The amplitude of the inhomogeneities (i.e., the contrast between the 
highest and the lowest baryon number densities) is determined by the value Vi of the wall 
velocity at T ~ T/v and by its minimum value Vm during the phase equilibrium stage. 

Obtaining quantitative estimates of the phase transition from analytical approximations 
is a difficult task. On one hand, the dynamics involves integro-differential equations which 
complicate the analysis. On the other hand, the exponential dependance of the nucleation 
rate with temperature introduces a large amount of uncertainty in the estimations. Nev- 
ertheless, an analytical inspection will provide a qualitative picture of the inhomogeneities 



and will be useful to study their parametric dependance. We dedicate the rest of this section 
to such a study. For that, we use the analytical expressions derived in Ref. [l^. In the next 
sections we will check these results by numerically integrating the equations for the progress 
of the phase transition. 



A. Dynamics of the electroweak phase transition 

Bubbles begin to nucleate immediately after T becomes smaller than T^. However, for T 
close to Tc the nucleation rate is vanishingly small, so the number of bubbles is insignificant. 
The onset of nucleation can be defined as the moment at which there are enough bubbles so 
that they begin to feel the presence of each other. More precisely, we define the temperature 
T/v as that at which the mean separation between bubbles equals the distance travelled by 
sound waves (which carry latent heat) since time tc Shortly after that, the temperature 
reaches its minimum value T^, say at time tm, and then increases due to the release of latent 
heat. Since the nucleation rate r(T) reaches its maximum at the minimum temperature T^, 
and is extremely sensitive to the temperature^, it is sharply peaked at t = t^- 

If the temperature T/v is close enough to the thin wall approximation can be used to 
estimate the nucleation rate. In this case, the number density of bubbles is given by 



The difference T^. — can be roughly approximated by'^ — T/v, where T/v is given by 

^' "-'^'K, (16) 



Tc — T]v / IGna^ 
and 

= 4 log {2^Mj,/Tc) + log (SL^T./STra^) + 6 log [(T, - T/v) /TJ (17) 

is a dynamical factor related to bubble nucleation at T ^ T/v. This factor depends on several 
parameters but its value is dominated by the first term, so ~ 100. 

Notice that the number of bubbles depends on the friction rj, and thus on the wall velocity. 
This is because for large latent heat is released more quickly, so the plasma reheats faster 
and the nucleation rate turns-off sooner. The friction coefficient depends on the viscosity of 
the plasma and the profile of the bubble wall. In the case (j)m(T) ~ T, which is required for 
baryogenesis, this dependence can be factorized in the form [13| 

r] ~ fjTa, (18) 

where is a dimensionless constant that depends only on the particle content of the plasma. 
The value of this parameter determines the initial velocity Wj. Using the approximation 



^ For T close to Tc, the dependance of F on T is dominated by an exponential function of the form 
F~exp[-C/ (Tc-T)']. 

^ In this approximation we are missing a slight (logarithmic) dependance of Tc—T„i on the friction coefficient, 
as can be deduced from Eq. (53) of Ref. 0|. 



V{T) ~ V (T/v) ~ L (T/v — Tc) in Eq. ()14|). we can write the initial wall velocity in the form 



which makes apparent a factor of v^^ in Eq. (jn 

As the Universe is reheated, bubble growth slows down. When bubbles have transferred 
enough energy to the plasma, the temperature reaches a value very close to Tc and stops 
increasing. This happens if the latent heat L is at least equal to the difference 6p = p (Tc) — 
p (Ttv), where 

p(T) = 7rVT730 (20) 

is the energy density of the plasma. If, on the contrary, L < 6p, then the temperature does 
not get close to and the expansion does not slow-down significantly. 

For L > 6p, the critical temperature is reached when the fraction of volume occupied by 
bubbles is 

/i ^ Sp/L. (21) 
The average bubble radius at the end of the reheating stage is thus 

^y^', (22) 



47rnfc/3 / 

Notice that in Eq. ()22j) the radius Ri does not have a straight dependance on the value of 
the wall velocity during this stage, as one could expect. This is because a larger implies 
a shorter reheating time. In fact, the time 6ti can be estimated as 6ti ^ Ri/vi. However, 
the bubble number density rib does depend on the initial value of v^, so in the end Ri does 
depend on Vi and Sti does not. 

After reheating, bubbles can release latent heat only at the rate at which this energy 
is taken away by the expansion of the Universe. Therefore, the rate of bubble expansion 
is given by the condition Lf ~ 4piJ, where H = ^ySrcGp/S is the expansion rate (here, 
G is Newton's constant). Since the fraction of volume in the broken-symmetry phase is 
f (t) ~ 4:TTnbR{t) /3, we can write / ~ ATmbR^Vm-, where the mean radius R2 during this 
stage is given by ribR^ ~ 1- Therefore, the wall velocity during slow combustion is given by 



(23) 



One could naively think that the velocities (fT^ and are unrelated, since Vi depends 
explicitly on the friction rj, whereas Vm is determined by the Hubble rate H. However, using 
the analytical approximations ()15II17|) . the velocities in the two stages can be written as 



We see that in this approximation Vm is also proportional to t]^^, so the ratio Vi/vf does 
not depend on i]. This is a consequence of the fact that in Eq. (jSHj), the number density 



of bubbles is proportional to rf, whereas p and H are not directly related to rj. Indeed, we 
have seen that rii, oc v^'^, so Vm is in fact proportional to f j, the proportionality factor being 
determined by the dynamics of reheating. From Eqs. ()24|) we obtain 



^ ^ 4 {L/T^f' K 

Since K and are both ~ 100, if L and a were of order 1 in temperature units, then 
Vi/vm ^ 10. However, in the case of the electroweak phase transition we have in general 
L/T'^ < 1 and a /T^ <^ 1, so Vi/vm can be quite large. 

We remark that the initial velocity is determined by the viscosity of the plasma and does 
not depend on the dynamics of the phase transition, whereas the ratio fj/fm depends only 
on the dynamics of the phase transition, so Vi and Vi/vm are quite independent quantities. 

A potentially important source of variation of Vw is bubble coalescence. When bubbles 
occupy a fraction of volume / ~ 1/3, they percolate. This means that at this point most 
bubbles are in contact, so they can group to minimize surface energy. This process con- 
tributes to the growth of bubbles and could dominate the wall velocity [l(. When / ~ 1/2 
coalescence stops because the regions of symmetric phase begin to form isolated bubbles sur- 
rounded by the broken-symmetry phase, so the interfaces are pushed again by the pressure 
difference between the two phases. Hence, there may be a different velocity variation for 
1/3 ^ / ^ 1/2. To ascertain whether this is the case we must compare the bubble growth 
rates caused by the two mechanisms. The growth rate due to coalescence is /c ~ {arii,/ p)^^"^ 
(l3l| . We compare it with the rate due to pressure difference during the slow combustion 
stage, /2 ~ pH/L. Therefore we have 



h ' ^T^gl'\ry^^) <'°"aT3 (^6) 

for Vi > 10^^. For values of L/T^ ~ 0.1 — 1, it would be necessary that a/T^ <^ 10^^ 
for coalescence to dominate and produce a significant departure from the velocity behavior 
described above. This does not seem likely for typical values of the electroweak phase 
transition parameters, although in general cr/T^ <^ 1^ 



B. Generation of inhomogeneities 

The amplitude of the inhomogeneities is determined by the total variation of the wall 
velocity during the transition. It is clear that larger density contrasts will arise if the 
velocity variation occurs far from the baryogenesis peak in Fig. Q i.e., in the region where 
the baryon asymmetry is most sensitive to f^. If Vi and Vm are much larger than fpeak, we 
can use the approximation ri^ oc v~^. If, on the contrary, -C fpeak, then the baryon 

density has a dependance oc f^. In any of these two cases, the ratio of the baryon 
densities in the high- and low- density regions, e = n^^^/n^", is given by the ratio vi/vm- 



^ In sectionl^we consider values of a/T^ in the range 10^'^ — 10^^. Notice that the values of L and a are 
not unrelated since both depend on the parameters of the theory. As a consequence, lower values of a 
correspond in general to lower values of L jl^] . 



In the intermediate case in which the variation of crosses the peak of n^, e may be much 
smaller. 

In general the peak velocity is in the range 10"^ < fpeak ^ 10~^ [lO, EH- Hence, to give 
the largest possible amplitude, the wall velocity should be either < 10~^ (in order to be 
^ 10~^) or Vw > 10~^ (so that Vw ^ 10~^). In general, however, the electroweak bubble 
walls have an initial velocity Vi ~ 10^^ — 10~^ 0, 113 5 so the velocity variation will not 
occur too far from the peak oi ub- If the initial velocity is ~ 10~^, then v^, will most likely 
cross fpeak duriug the transition. So, the most favorable situation is that in which Vi is closer 
to 10~^. In this case the velocity will not be far from Vpeak initially, but will probably not 
traverse it either. Large inhomogeneities can then arise if the minimum velocity departs 
enough from Vi. 

As we have seen, the velocity variation depends on L and a. According to Eq. (j25|) . for 
large L or small cr, the wall velocity will change considerably, and important baryon density 
contrasts may be formed. On the contrary, if L is too small or a too large, the wall velocity 
will not change significantly during the transition. In this case, the baryon number density 
will be rather homogeneous. Notice that for L — >■ there is no reheating at all, and we 
should have Vi/vm 1, so Eq. ()25p fails in this limit. The problem is that for very small 
L, the bubble expansion never enters the slow combustion stage, so the basic assumptions 
that lead to Eq. ()23p for the minimum velocity break down. 

Since all the bubbles nucleate in a short time 5tr at the beginning of the transition, all 
the inhomogeneities should have the same size scale, determined by the separation between 
centers of nucleation, d . As can be seen in Eqs. fll5lll8|l . rif, is a very sensitive 

function of L, a, and the temperature difference Tc — T^,, which in turn depends on L and a 
too. Furthermore there is an additional dependance on a through the friction coefficient rj. 
So, we cannot trust the dependance of c? on L and a obtained from these approximations. 
On the other hand, we do not expect a significant dependance of on the friction, so we 
can rely on the behavior rib oc f which implies d oc Vi. 

The spherically symmetric inhomogeneities that originate inside the bubbles can have 
two basic geometries, depending on the relation between fpeak, Vi, and Vm- If fi < fpeak, the 
baryon asymmetry generated at the bubble walls decreases with time as bubbles expand. 
Hence, a larger baryon number density is produced at the centers of bubbles, and then the 
geometry of the inhomogeneities is that of spheres of high baryon number density ~ (fj), 
surrounded by regions with lower density ~ (f^) (see Fig. These spheres have radius 
Ri and occupy a fraction of volume /i, given by Eqs. (pij) and (HH). 

In the case Vm > fpeak, the baryon production will be larger at the end of bubble ex- 
pansion, and inhomogeneities will have the form of walls which surround spherical voids of 

— 1/3 

radius Ri. Therefore, these walls have a width w ^ {fi — l)Ri and occupy a volume 
fraction /2 = 1 — /i- Of course, there is a third possibility, namely, that Vi > fpeak and 
fm < fpeak, which case the structure of the inhomogeneities can be more complex. How- 
ever, as explained above, we do not expect inhomogeneities with large amplitude when f^, 
varies across fpeak, so this case is in fact irrelevant. 

According to Eq. (PT|) . the volume fraction fi is given by the ratio Sp/L, which must be 
less than 1 in order to have acceptable reheating and velocity variation. By naturalness, L is 
not likely to be too close to Sp, so /i is not expected to be very close to 1. This implies that 
the width w should be at least of the same order of Ri. On the other hand, /i could be very 
small if L ^ 6p. In this case, we would have inhomogeneities of a small size Ri ~ (fi/nt,)^^^, 

— 1/3 

separated by wide walls of size w ^ d ^ . Notice however that Ri <ti d is not hkely. 





FIG. 2: Schematic picture of the geometry of the baryon inhomogeneities created during bubble 
expansion. Shaded regions correspond to high baryon number density. The left figure shows the 
case in which Vi < Vpeakj whereas the right figure shows the case Vm > t'peak- 

Indeed, a size Ri ~ O.ld already requires a latent heat three orders of magnitude larger than 
the energy difference between T/v and Tc. Although 6p/L can change significantly from one 
model to another we do not expect that this ratio can be made arbitrarily small in a 
realistic theory, since both the latent heat and the amount of supercooling depend on the 
effective potential. As a consequence, we expect that the sizes of the spheres and walls in 
Fig. 121 will be similar for reasonable values of the parameters. This is important because if 
the size of the inhomogeneities is too small, they will be quickly erased by diffusion and will 
not survive until later epochs. 



IV. NUMERICAL CALCULATION OF THE THE ELECTROWEAK PHASE 

TRANSITION 

In order to simplify our numerical calculation of the phase transition, we will use a free 
energy density of the form 

V {(f), T)=D (T^ - T^) <f - ETcf + ^<j)\ (27) 

It is well known that to one-loop order in perturbation theory and in the high-temperature 
approximation, the effective potential takes this form, where the parameters D, Tq, E, and 
A depend on the particle masses (see e.g. jlH). The actual form of the free energy depends 
on the model and differs in general from ()27|1 . This may happen, for instance, in the two- 
loop approximation, or in the case perturbation theory or the high-temperature expansion 
are not valid. In any case, this approximation is useful to simulate the phase transition in 
different models, which is convenient for a general treatment. 

This is accomplished by adequately choosing the parameters in Eq. (P7j) so that the free 
energy gives the correct values for the thermodynamic parameters defined in the previous 
section. For instance, the critical temperature and latent heat are given by 

T, = To/Vl - E^XD, 

L = 8D{E/XfT^Tl (28) 

The first-order phase transition can take place between the temperatures Tc and Tq, since 
in this temperature range the potential ()27p has a local minimum at = 0, corresponding 



to the symmetric phase, and a global minimum at 



0m (T) 
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(29) 



corresponding to the broken-symmetry phase. These two minima are separated by a barrier. 
At T = Tc the minima are degenerate, and at T = Tq the barrier disappears and the 
point = becomes a maximum. The pressure difference between the two phases is 
V (T) = V (0m (T) , T). The order parameter of the phase transition is given by 



0,/T, = 0„ (Te) /Te = 2E/X. 



(30) 



The progress of the transition is characterized by the fraction of volume occupied by 
bubbles of the broken-symmetry phase. In the case of the electroweak phase transition, this 
is given by 0, H 13 



fit) 



1 - exp I -y ^* r (T') R [t', dt'l . 



(31) 



Here, R {f, t) is the radius of a bubble that nucleated at time t' (and temperature T') and 
expanded until time t, 



R{t',t) 



{T") dt" 



(32) 



where the wall velocity is given by Eq. (fT^ . and we have neg lected the initial radius of the 
critical bubble. The nucleation rate F (T) is given by 



(33) 



where (T) is the three-dimensional instanton action^, which coincides with the free energy 
of a critical bubble in unstable equilibrium between expansion and contraction. 



An I r'^dr 
'o 



(34) 



The configuration of the nucleated bubble can be obtained by extremizing this action. 
Hence it obeys the equation 

dU 2d(f) dV 



dr"^ 



r dr 



For T ^ Tc the radius of the bubble becomes infinite, so the profile of the wall can be easily 
calculated by neglecting the second term in Then, we can readily compute the wall 

tension for the model ^T7\\. 

^ 2v^ 3 

^ 3A5/2 



(36) 



^ The prefactor of the exponential in is roughly assumed to be of order T^, since the rate is dominated 
by the exponential. 



This equation, together with Eqs. and (jHUj) determine completely the values of the 

parameters of V (0, T) in terms of thermodynamical properties of the free energy. 

Analytical approximations for the profile of the critical bubble introduce large errors in 
the nucleation rate. Therefore, we will use the numerical fit for 5*3 given in Ref. (27l |. 



f = 13.72^ a'/'g (a) , (37) 

, a A 2.4 0.26 

4 V I - a {I - a) 

where a = XD (T^ — Tq) /E^T"^. This dimensionless parameter can also be written as 



2 



where it is apparent that it varies from 1 to as T decreases from Tc to Tq. 

Finally, the variation of temperature with time during the phase transition is given by 



^ " a3 +27rV/45^- ^^^^ 

Here, the prime means total derivative with respect to temperature. The first term in 
Eq (j39|) gives the usual variation of T in the adiabatic expansion, while the second term 
accounts for the release of entropy during the phase transition. The evolution of the scale 
factor a{t) is given by the Friedman equation, [a/ of = H"^ = 87rGp(T)/3. Since the 
duration of the electroweak phase transition is much shorter than the age of the Universe, 
we can use the approximation if ~ if^ [13 ■ Note that in Eqs. (jHTj) and (jH^ we have 
neglected the variation of length scales due to the expansion of the Universe. We cannot do 
the same in Eq. (jHIJjl . since for T ~ small changes of temperature are relevant. Indeed, 
the scale of temperature variation in a first order electroweak phase transition is given by 
T,-To~(EV2AD)T,<T,. 



V. NUMERICAL RESULTS 

In order to solve numerically Eqs. ()31II39|1 it is convenient to differentiate Eq. (j39j) . At 
this stage we can make the approximations a ~ Oc, T ~ Tc, so we have 



15 V" 
2^r2 



/ 



dT _ Idt 15 V 
7^ ~ ~2 tT ^ 2^T3^' 



In addition, the variation of / is given by 



(40) 



df = il-f)A7rv^ (t) 



TR^dt' 



dt. 



(41) 



Integration of these equations gives the values of the radius R{t) and wall velocity Vw{t) at 
every time. Then, the profile ^^(r) is obtained with the aid of Eq. (fTUj) . as explained at 
the end of section |n] For the numerical calculation we have made a simple discretization 
of Eqs. ()40|) and ()41|) . We have used a time step much shorter than the duration of the 



phase transition (~ 1/1000), and we have verified that the result remained unchanged as 
this step decreased further. Therefore, the uncertainties of the numerical calculation are 
certainly smaller than the ones introduced by the analytical approximations used in the 
previous section. 



A. Parameter values 

Several extensions of the Standard Model provide a strong phase transition as well as 
enough CP violation for baryogenesis. Supersymmetric models were extensively studied in 
the literature, especially the case of the MSSM, which gives the required BAU in a light 



right-handed stop scenario [10, IJJJ, |2J, |2a |22|- Some alternatives to this model include 
non- minimal supersymmetric extensions (see e.g. ^^), the addition of heavy fermions to 
the SM js^, a dimension-six operator in the effective Hi ggs potential |4o| . or the presence 
of hypermagnetic fields during the phase transition |4lll4^. Therefore, we will consider 
various values of the parameters. 

For the case of the MSSM, the thermodynamic parameters of the phase transition have 
been calculated Hence, we will take this starting point for the parameter 

variation. We will consider different values of L and a, since these two parameters are the 
ones which most directly affect the dynamics of the phase transition. Thus, according to the 
non-perturbative study of Ref. our reference values will be Tc ~ 85GeV, (pc/Tc — 1, 

L/T^ ~ 0.4, <j/T^ = 0.01. The viscosity of the plasma has also been studied for the MSSM 
in the light stop scenario i^]. We will consider values of the friction parameter fj ~ 0.1 — 10, 



which give initial velocities^ Vi ~ 10~^ — 10~^, in accordance with the results of Ref. [2^ 

To illustrate the kind of profiles of the baryon number density that can be created inside an 
expanding bubble, we will consider several values of the electroweak baryogenesis parameters. 
According to the discussions in section llllj the profile depends essentially on the position 
of fpeak relative to the velocity variation. We have seen in section |n] that setting the values 
of the parameters a and 6 to ~ 1, the peak velocity is given by fpeak ~ y/DT^s- For the 
MSSM, D ~ 100T~^, so this gives the correct value Vpeak ~ 10~^. In other extensions of 
the SM the effective diffusion constant D may differ from that value, and further deviations 
may arise through the precise values of a and b. Therefore, we will consider values of fpeak 
in the range 10~^ — 10~^. 



B. Dynamics of the phase transition 

The evolution of the phase transition for these values of the parameters is shown in Fig. 
El We have plotted the variable a defined in Eq. flHHI) . the fraction of volume / occupied 
by bubbles, and the bubble wall velocity v^, as functions of time. Notice that the evolution 
of the temperature and the fraction of volume is affected by the bubble wall friction only 
at the beginning of the phase transition. This is because during the reheating stage, the 
bubble expansion is dominated by friction. On the contrary, during the phase equilibrium 
stage the development of the phase transition is determined by the balance between the rate 



^ We obtained these relations between the values of ij and Vi numerically. The values of Vi obtained from 
Eqs. (|24|l and l|18|l are a factor ~ 2.5 less than the numerical values. 



of latent heat release and that of cooling due to the expansion of the Universe, as explained 
in section ITTTl 
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FIG. 3: Fraction of volume, temperature parameter a, and wall velocity as functions of dimension- 
less time T = {t — tc) /tc- 

We see that the wall velocity decreases by a factor that does not depend on the friction 
coefficient (and hence on the initial velocity Vi), as anticipated by the analytical approxima- 
tions. In the present case the factor is ~ 30, to be contrasted with the estimation (j2Sl), which 
for the parameters under consideration gives Vi/vm — 100. As expected, the discrepancy is 
of order 1. 

Notice also that the wall velocity is not constant in the slow combustion stage, but it 
slowly begins to grow again after reaching its minimum value Vm- This happens because as 
the regions occupied by the symmetric phase become smaller, the energy must be released 
more quickly at the interfaces in order to compensate the cooling produced by the adiabatic 
expansion. 

In the analysis of the inhomogeneities we have assumed that all bubbles begin to expand 
at the time tj at which the nucleation rate F takes its maximum value. We have checked 
numerically that most of the bubbles nucleate at the very beginning, while the volume 
fraction increases from 10~^ to 10~^. Moreover, Fig. |3] shows that F is active during a time 
interval (5tr which is two orders of magnitude shorter than the total duration of the phase 
transition (compare with Fig. This introduces a dispersion of about a 5% in the bubble 
radius at the end of the transition, which we define to be the moment at which / = 0.99 
(the initial and final times delimit the plot of the wall velocity in Fig. Oj). Hence, there will 
not be inhomogeneities at different scales, although some dispersion in sizes will arise due 
to the fact that the centers of nucleation are randomly scattered throughout space. 

C. Profile of the baryon inhomogeneities 

In Figs. El and IHl we show the inhomogeneity profiles that arise for initial wall velocities 
slightly above 0.1 and slightly below 0.01 respectively, and different possible values of fpeak- 
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FIG. 4: The nucleation rate as a function of dimensionless time r defined in Fig. |S1 
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FIG. 5: Baryon number density as a function of the bubble radius r in units of the age of the 
Universe tc, for Vi ~ 0.1. 

In the case Vi > 0.1, we see that larger amplitudes are obtained for smaller values of 
■^peak, i-e., when Vi is more distant from the baryon generation maximum. The profile for 
^^peak ~ 10^^ in Fig. El is in good agreement with the result of Ref. where a dependance 
hb oc was assumed. The maximum of this curve corresponds to the minimum wall 
velocity Vm reached during the phase transition, since in this case the velocity fpeak is not 
attained. In the other curves, instead, the peak in ^^(w^) becomes manifest. In the case 
Vpeak ~ 10^^, the amplitude is smaller, and we observe the presence of two hills in the 
curve. This is because crosses Vpeak immediately before reaching the minimum Vm- Then 
it crosses the peak again due to the final acceleration. The case fpcak ~ 10~^ gives still a 
smaller amplitude, since the velocity variation is more symmetric around the peak. 

Fig. M shows the case Vi < 10^^. As expected from the discussion of section IIIH the 
profiles behave quite oppositely to the previous case. The largest amplitude is obtained for 
Wpeak ~ 10~^, and in this case the baryon number density is higher at the centers of the 
bubbles. The upper curve is similar to the lower curve of Fig. El because in this case the 
wall velocity crosses the peak. We see that when this happens, the amplitude is smaller and 
the baryon inhomogeneities are not clearly localized either in the interior or the exterior of 



the spheres of Fig |21 
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FIG. 6: Baryon number density as a function of the bubble radius r in units of the age of the 
Universe tc: for Vi ~ 0.01. 



It is important to notice that these profiles describe the inhomogeneities produced by 
an average bubble that expands without colliding with other bubbles. Furthermore, bubble 
collisions may have some effect on the development of the phase transition. Since Eq. ()31|) 
takes into account overlapping of bubbles, the present simulation of the phase transition is 
in principle reliable in the whole time interval, with the possible exception of collisions and 
coalescence, which occur when the fraction of volume is / ~ 0.3 — 0.5. We have seen in 
section IIIII that the dynamics of the electroweak phase transition is not affected by bubble 
coalescence. However, we are neglecting some possible effects, such as additional reheating 
due to bubble collisions. Such effects could introduce some distortion in the profiles shown 
in Figs. Elandini 



D. Size and amplitude of the inhomogeneities 

The size scale of the inhomogeneities is given by the distance between centers of bubbles, 
d ~ n^^^^ We plot this distance in Fig. [3 as a function of the bubble wall tension for 
different values of the latent heat and friction. Notice that, as expected, changing r/ by an 
order of magnitude induces the same variation in d. We have studied the dependance of 
d on the initial wall velocity in the range 10~^ < Vi < 0.4, and we have checked that the 
linearity d oc is verified for the different values of L and a. 

The maximum possible baryon density contrast e is determined by the ratio Vi/vm- In 
Fig. IHlwe plot this ratio as a function of latent heat, for different values of a. For MSSM 
parameters we find Vi/vm ~ 10 — 100. We observe that for large L, Eq. ipHj) gives the correct 
order of magnitude and the exact parametric dependance. In contrast, for smaller values of 
latent heat, the behavior of the numerical curves changes abruptly. This was expected, since 
when L < 6p the reheating becomes insufficient and the phase transition ceases to have two 
well defined stages. Thus, the curves depart from the approximate behavior ()25|) . and for 
L ^ 0, Vjn approaches required physically. 
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FIG. 7: Size of the inhomogeneities as a function of the bubble "wah tension for different values of 
the latent heat and friction. 
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FIG. 8: The velocity ratio Vi/vm as a function of latent heat for different values of the bubble wall 
tension. 

In this "work we have assumed that the latent heat released as the bubbles expand is 
quickly distributed, so that the only effect is a global reheating. To estimate whether 
this approximation is correct or not, we must compare the time scales for the evolution of 
the phase transition with the time in which the latent heat liberated in a deflagration front 
reaches the neighboring bubbles. As seen in Fig. 01 for our reference values of the parameters 
the quick reheating stage occurs in a time 6ti ~ 10~Hc. On the other hand, the time it takes 
a shock front to travel the distance between bubbles is given by 6tsh ~ d. Hence, we expect 
that our treatment will be reliable for d/tc <^ 10~^. According to Fig. [3 for Vi > 10~^ 
{fj = 0.5) this definitely requires a/T^ < 10~^. For smaller values of the initial wall velocity 
this condition is relaxed. After the system reheats, the transition proceeds more slowly. In 
this stage the time scale is 6t2 ~ 10~^tc; and the approximation breaks down only for higher 
values of the wall velocity and tension, Vi ~ 0.1, a/T^ ~ 10~^. Since the amplitude of the 
inhomogeneities is larger for smaller values of a (see Fig. IH)), in the interesting cases the 
global reheating approximation is valid. 



VI. EFFECT ON THE QUARK-HADRON PHASE TRANSITION 



The aim of this paper was to study the characteristics of the baryon inhomogeneities 
generated in the electroweak phase transition. However, we wish to comment briefly on the 
possibihties for these inhomogeneities to leave a sequel in the subsequent evolution of the 
Universe. The time scale for the wash-out of the inhomogeneities depends on their size and 
amplitude, and on the diffusive processes. For temperatures in the range lOOGeV^ < T < 
IMeV, the only important dissipation process is neutrino diffusion. A wide range of sizes 
and amplitudes is left almost unaffected by this process i3- Therefore, it is likely that 
electroweak inhomogeneities survive somewhat unchanged until the QCD scale 0, 17 1. in 
which case they may influence the dynamics of the quark-hadron phase transition. Below 
IMeV other processes (namely, baryon and photon diffusion) become important and may 
completely erase the electroweak inhomogeneities before the nucleosynthesis epoch. 

The evolution of baryon inhomogeneities between the epochs of T = lOOGeV and 
T = lOOMeV has been considered in Ref. According to this calculation, a baryon 

inhomogeneity which at the electroweak scale has amplitude < 10^ and size > 10~^tc, corre- 
sponding to the cases considered in Figs. [7| and |Hl may survive almost unchanged until the 
QCD epoch. 

As shown in Fig. |S1 in the MSSM inhomogeneities with amplitude enssM ~ 10 — 100 are 
like ly. However, there exist several interesting scenarios for the electroweak phase transition 
MMWd, 41, 4J besides the MSSM. In some of them the phase transition is much stronger 
than in the latter. In such extensions of the SM, the parameters affecting the dynamics will 
differ significantly from those of the MSSM. In particular, larger values of the latent heat 
are expected in stronger phase transitions This leads in general to a larger e, as can 
be seen in Fig. |H1 Further enhancement (or suppression) may also result from the variation 
of the wall tension. In fact, one expects larger values of a in stronger phase transitions, so 
in principle this effect goes in the opposite direction to that of latent heat. In any case, a 
sensible estimation of the inhomogeneity amplitude requires a calculation of L and a in each 
particular case. Nevertheless, it is reasonable to expect that in some models amplitudes 
considerably larger than in the MSSM may arise. In such models we will possibly find 
e > 100. 

As can be seen in Figs. and |H1 values of L and a which give larger amplitudes, yield 
also smaller sizes. Such a correlation is important because inhomogeneities with larger 
amplitudes and smaller sizes are more easily diluted. Extrapolating our results beyond the 
ranges of parameters considered in the previous section, we see that for an amplitude ~ 10^ 
we would have a size d/tc ~ 10~^ — 10~^, which is still undamped by neutrino inflation. 
Larger amplitudes will probably be affected. For instance, according to the results of Ref. 
[l7l |. an initial amplitude ~ 10^ may decrease to ~ 10^ if the size scale is < 10~^tc- Therefore, 
it seems unlikely that inhomogeneities generated in the electroweak phase transition may 
have an am plit ude much larger than 10^ at the time of the QCD phase transition. 

In Ref. it was found that a significant effect on the latter is obtained for inhomo- 
geneities with an amplitude e ~ 10^. However, it was also noticed that smaller values of 
e will have similar effects if the amount of supercooling is smaller than the one considered 



there, in accordance to Ref. |4^. Furthermore, in the treatment of Re f. i lQ] the possibility 
that the nucleation rate changes with the chemical potential (see e.g. j45|) was not consid- 
ered. Taking into account this effect could affect strongly the discussion in Ref. and 



therefore lead to a significantly different value of the inhomogeneity amplitude needed to 



affect tliis transition. Altogetlier, it is clear tliat tlie dynamics of tlie quark-liadron pliase 
transition is still not fully understood, and we can sensibly expect that it may be affected 
by prior inhomogeneities with amplitudes of the order of those that are generated in the 
electroweak phase transition. We hope to address this issue in future work. 



VII. CONCLUSIONS 

In this paper we have made a detailed study of the baryon number inhomogeneities 
that can be generated as a byproduct of electroweak baryogenesis. In particular, we have 
investigated analytically the parametric dependance of the inhomogeneities, and we have 
made a more precise calculation of the amplitude and size of the density contrasts by solving 
numerically the equations for the development of the electroweak phase transition. By 
doing that, we have checked several analytical approximations derived in Ref. 0] for the 
dynamics of the phase transition. An important difference between our treatment and 
previous analysis is that we have taken into account the fact that in general the baryon 
density has a maximum at a bubble wall velocity = fpeak- We have accomplished this by 
deriving a simple analytical approximation for the dependance of the baryon asymmetry of 
the Universe on Vw 

The characteristics of the inhomogeneities depend on the parameters that determine 
the dynamics of the phase transition, namely, the latent heat L, bubble wall tension cr, 
and friction coefficient 77, and on the parameters that govern the production of baryons, 
which are essentially the diffusion constant D for the chiral quark asymmetry and the weak 
sphaleron rate T^g. Some general features, however, are independent of these parameters. 
For instance, the spherical symmetry and the absence of a significant dispersion in sizes. 

Indeed, we have shown that all the bubbles are nucleated in a short time interval at 
the beginning of the transition, and therefore have roughly the same radius. (We wish to 
remark here that this feature of first-order phase transitions may be important for other 
cosmological consequences as well, such as the formation of topological defects or magnetic 
fields. We will address these issues somewhere else.) Since the baryon inhomogeneities are 
generated in the walls of expanding bubbles and depend on the wall velocity, this implies 
that all inhomogeneities have approximately the same size, amplitude, and profile. 

The size scale is given by the distance between centers of nucleation, so it depends only 
on the phase transition parameters L, a, and rj. In particular, we have seen that the higher 
the initial wall velocity, the larger the bubble separation, due to the sooner turn-off of 
the nucleation rate. The amplitude e of the density contrasts may depend further on the 
baryogenesis parameters; more precisely, on Vpeak — \fDTws- However, if Vw stays far from 
fpeak throughout the phase transition, e is given by the ratio Vi/vm between the maximum 
and minimum values of Vw We have seen that this ratio does not depend on the friction 
parameter, so in this case e is only determined by L and a. If the wall velocity crosses the 
value fpeak during the transition, then the amplitude is smaller than Vi/vm- 

The exact profile of the inhomogeneities, on the other hand, depends both on Upeak and 
on the dynamics of the phase transition. In the case of the MSSM, the initial velocity is 
most hkely close to the baryogenesis peak, Vi ~ fpeak ~ 10~^ [H, l^l- As we have seen, this 



implies that the geometry is that of high-density spheres surrounded by low-density walls. 
This is contrary to the case of the SM, where wall-shaped inhomogeneities are formed j^, 
showing that both kinds of profile can arise in different models. 

For MSSM values of the parameters L and a, we find amplitudes e < 100. Nevertheless, 



since the ratio Vi/vm grows with L^/^ and with a"^, e may be much larger in other extensions 
of the SM. We have seen that if this is the case, the electroweak baryon inhomogeneities will 
probably affect the dynamics of the quark-hadron phase transition, as discussed in Ref . 

For the size of the inhomogeneities, we found in the case of the MSSM a value of order 
10^^—10"^ times the age of the Universe (in agreement with Ref. However, this quantity 
is sensitive to r] and a, and can deviate significantly from this value in other models. 

These behaviors are illustrated in Figs. [7| and |H1 Our analytical expressions agree with 
the numerical computations within the order of magnitude. We have considered parameters 
that varied in a wide range in order to cover possible extensions of the SM, other than the 
MSSM. Furthermore, our results can be easily extrapolated beyond these ranges. 



Acknowledgements 

A. M. is grateful to M. Quiros for useful discussions, and to the Physics and Mathematics 
Institute at Michoacana University, where this work was initiated, for kind hospitality. The 
work by F. A. was supported by CIC-UMSNH and Conacyt grant 32399-E. 



APPENDIX A: THE CHIRAL ASYMMETRY 

In this appendix we find an analytical approximation for the density nL{z). In section 2 
we have seen that this is of the form 

riL = Ae-'^»+^ (Al) 

and we argued that A (x v^, which is all we need to know for the present paper. However, this 
may not be true in general, so it is important to see under which conditions the coefficient 
A depends only linearly on t>^. 

As we have seen, the solution to Eq. (jS)) is given by Eqs. (jlj and with different values 
of r, and consequently of uj±, A, and B in the broken and symmetric phases. Specifically, 
Ff, = F, Fs = aTy^s, where the index b (s) stands for the broken (symmetric) phase. According 

to Eq. (0) we then have uji± = Vw/2D ± [(t>^/2D)^ + Ti/ ^, and the constants Ai and 
Bi are given by the boundary conditions as explained in section |Tl| As a result, one obtains 

/■O 

A, = [D {uj,+ - u,^)]'^ / €^"+^7 {z) dz, 

J —oo 

POO 

Bs = [D {ujs+ - ujs_)r' / (^) dz, (A2) 

Jo 

and 

_ (^b+ — ^ ^s- — ^b- p 

— ^b -Ds; 

"-^8+ ~ ^b- ^^s+ — ^^b- 

— ^b- ~ ^b- 

The source 7 is proportional to the wall velocity, 7 = ^^7. If 7 {z) is localized in some 
region around the wall, then it can be easily seen that the coefficients flA2|) are such that 



outside this region only the exponential with the right sign of uo survives in each phase. To 
obtain an analytical result, we need to use an approximation for 7. The simplest approxi- 
mation for a function 7 [z) that is localized inside the wall is a step function 2J|, say, 7 = 70 



for —Lyj < 2; < 0, and 7 = outside. (The value of the constant 70 depends on the CP 
violating force inside the wall.) Then, the solution for z > is of the form (jAljl . with 

A = '^-7"''>"-^°. (A4) 

Notice that for f ~ lO^^T, the relation (f^/2D)^ ^ Tb/D holds for any possible value 
of the wall velocity, as long as D > lOT^^. So, the u!b± are insensitive to the wall velocity. 



LUb± — ±yT/D. Furthermore, since Tg <C Tb, we can also neglect uig^ in Eq. ()A4jl . and the 
coefficient A becomes 

A ~ ^ ^, (A5) 

which is evidently linear in Vy^. This expression for the coefficient A coincides with the one 
obtained in Ref. ji^. On the other hand, the 2r-dependance in Eq. (jAlll is different due to 
our inclusion of the sphaleron process in the diffusion equation for rii- As we have seen, this 
modification does not introduce any qualitatively different behavior in the final result for 
the baryon number density. However, some quantitative 0{1) difference in the BAU shows 
up for vl < DT^s, i.e., for < 10"^ 
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